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1.  INTRODUCTION 


Electromagnetic  fields  incident  on  a  scattering  object  induce  electric  and/or 
magnetic  currents  which  radiate  to  produce  scattered  fields.  An  image  is  a  syn¬ 
thetic  picture  of  an  object  generated  from  electromagnetic  scattering  data  rather 
than  the  direct  observation  of  the  target  [1],  The  IEEE  definition  for  image  is  “a 
spatial  distribution  of  a  physical  property  such  as  radiation,  electric  charge,  con¬ 
ductivity  or  reflectivity,  mapped  from  another  distribution  of  either  the  same  or 
another  property”  [2].  An  image  permits  human  visual  recognition  and  identifica¬ 
tion  of  the  object  which  may  also  carry  detailed  information  about  the  scattering 
characteristics  of  the  object.  Imaging  techniques  have  a  large  variety  of  applications 
such  as  target  recognition  and  identification,  diagnostic  methods,  acoustic  imaging 
for  medical  diagnosis,  seismic  imaging  and  many  others. 

Most  of  the  image  reconstruction  techniques  are  based  on  the  fact  that  a  two-di¬ 
mensional  (2-D)  image  of  an  object,  which  can  be  represented  by  a  2-D  distribution 
of  scattering  centers,  and  its  2-D  frequency  spectra  are  Fourier  transform  pairs 
[3,4,5] .  Therefore  the  image  of  a  2-D  object  can  simply  be  generated  by  inverse 
Fourier  transforming  the  2-D  frequency  domain  data.  The  imaging  algorithm  used 
in  this  paper  utilizes  1-D  filtered  time  domain  signals  together  with  aspect  angle 
information  and  replaces  the  conventional  2-D  Fourier  transform. 

The  scattering  characteristics  of  most  targets  are  polarization  dependent  which 
makes  the  images  polarization  dependent  also.  Hence,  polarization  is  an  impor¬ 
tant  parameter  for  target  identification  by  imaging.  In  this  report,  images  of  both 
geometrically  simple  (like  sphere,  ellipsoid,  square  plate)  and  complicated  (like  an 
aircraft)  targets  using  different  polarizations  are  generated  and  compared.  Two 
different  sets  of  colors  assigned  to  vertically  (VP)  and  horizontally  (HP)  polarized 
images  in  order  to  display  the  polarization  dependence  of  the  images. 
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Scattering  measurements  of  several  targets  were  made  using  The  Ohio  State 
University  Compact  Radar  Cross  Section  Measurement  Range  [6].  Most  of  the 
measured  data  span  the  range  2-18  GHz.  For  some  of  the  targets,  the  frequency 
scaling  technique  permitted  data  between  0.01-18  GHz  to  be  obtained.  Calculated 
data  were  generated  using  approximations  based  on  uniform  asymptotic  techniques 
in  terms  of  the  Geometrical  Theory  of  Diffraction  (GTD)  [7].  The  same  frequency 
range  was  used  for  measured  and  calculated  data  in  order  to  directly  compare 
measured  and  calculated  images. 

Chapter  2  covers  the  basic  theory  related  to  the  image  reconstruction.  Relation¬ 
ship  between  the  frequency  spectra  and  the  image,  imaging  algorithm  used  and  the 
sampling  criteria  are  all  discussed  in  this  chapter.  Chapter  3  covers  the  applications 
of  the  theory  discussed  in  Chapter  2.  Also,  measurement  and  data  processing  tech¬ 
niques  used  are  discussed  and  some  results  using  either  measured  or  simulated  data 
with  different  polarizations  are  illustrated.  Conclusions  drawn  from  the  research 
are  discussed  in  Chapter  4. 

2.  BASIC  THEORY 

The  image  of  an  object  can  be  characterized  by  a  2-D  distribution  of  the  scat¬ 
tering  centers  which  is  related  to  the  2-D  frequency  spectra  by  a  2-D  inverse  Fourier 
transform  relationship.  This  relationship  has  been  derived  and  discussed  in  Section 
2.1.  Section  2.2  covers  the  imaging  algorithm.  Limitations  on  the  images  due  to 
space  and  frequency  sampling  criterion  are  discussed  in  section  2.3. 


Figure  2.1:  Imaging  geometry  for  a  two-dimensional  object. 


2.1  Image  Reconstruction  Using  2-D  Frequency  Spectra 
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To  understand  how  an  image  is  generated  using  either  measured  or  simulated 
data,  assume  a  rotating  object  which  is  uniformly  illuminated  by  a  stationary  trans¬ 
mitter  as  shown  in  Figure  2.1.  The  object  is  rotated  through  an  angle  8  and  scat¬ 
tered  signals  are  recorded  as  a  function  of  frequency  and  aspect  angle.  Also  assume 
that  the  source  is  monochromatic  (CW)  with  a  frequency  /.  The  objective  is  to 
obtain  a  2-D  image  of  this  object  in  the  cross-sectional  plane.  The  cross-section  of 
the  object  is  represented  by  a  2-D  distribution  of  scattering  centers,  g(x,y),  which 
is  denoted  as  the  “reflectivity  density  function”  [4,5].  Assume  y  and  x  are  the  fixed 
and  u  and  v  are  the  rotated  down  and  cross  range  coordinates,  respectively.  For  a 
fixed  8 ,  the  total  received  signal  will  be  the  summation  of  the  contributions  of  all 
the  scattering  centers;  i.e.,  the  integral  of  g[x,y)  over  u  and  r  assuming  that  for 
a  fixed  range,  all  scattered  fields  are  in  phase.  A  phase  factor  of  e~^7kv  must  be 
supplied  to  account  for  the  phase  delay  along  the  down-range. 


(2.1) 


G(B)  =  r  r  ge{u,v)e-]2kvdu  dv 

J  —  oc  J  —  oo 

=  f  [  gt{u,v)e~j4*vlxdu  dv , 

J  —  oo  J  -oo 


where  k  is  the  wave  number  and  50(1/, v)  is  the  reflectivity  density  function  for  a 
particular  aspect  9. 

Rotated  and  fixed  coordinates  are  related  as, 


cos  9  —  sin  9  a 

sin  9  cos  9  v 


(2.2a) 


u  cos  9  sin  9  x 

v  —  sin  9  cos  9  y 

Using  Equation  (2.2b)  in  Equation  (2.1),  one  obtains 


/oo  yoo 

/  g{x, y) 

-00  J  -  00 


and  defining  the  new  variables  /*  and  fy  as, 

r  ^  2  si  n9 

J*  =  \  * 


(2.2b) 


(2.4a) 


a  2  cos  9 


(2.4b) 


Equation  (2.3)  becomes, 


£'(/*,/„)=  I"  r  9(*,y)r-j2'U'x*f'y)drdy, 

J  —  00  J  -00 


which  is  a  Fourier  transform  relationship.  Therefore,  g(x,y)  can  be  reconstructed 
by  inverse  Fourier  transforming  (*'(/*, /v), 

g(T,y)Ui-G(fx,fv),  (2.6) 


9ir,y)=  r  r  Gif,Jv)ei7'U'r*f"')dfrdfv. 

J  -  oc  j  -  00 


4 


With  previous  assumptions,  g(x,y)  of  Equation  (2.7)  corresponds  to  a  monochro¬ 
matic  wave  of  frequency  /.  For  a  band  of  frequencies,  the  resultant  <7(x,j/)  must  be 
the  summation  of  the  contributions  of  each  frequency  component. 

9(r,y)  =  (2.8) 

1=1 

or,  using  the  linearity  property  of  the  Fourier  transform, 

0(*,y)  =  F~'  |^G(/r.,/fc)|  ,  (2.9) 

which  is  a  2-D  Fourier  transform  relationship.  Therefore,  one  can  recover  g{x,y) 
simply  by  inverse  Fourier  transforming  the  2-D  frequency  spectra. 

2.2  Imaging  Algorithm 

In  Section  2.1  it  was  shown  that  the  image  of  an  object  and  its  2-D  frequency 
spectra  are  Fourier  transform  pairs.  Hence,  the  image  of  a  target  can  be  recon¬ 
structed  by  simply  inverse  Fourier  transforming  the  2-D  frequency  domain  data. 
On  the  other  hand,  measured  data  are  usually  sampled  on  a  polar  raster  rather 
than  a  cartesian  one  and  need  to  be  converted  to  a  cartesian  raster  in  order  to  get 
compatible  data  for  the  2-D  FFT  algorithms.  Conversion  from  polar  to  rectangular 
coordinates  requires  interpolation  of  the  data  points  which  makes  quality  of  the 
final  image  very  dependent  on  the  method  of  interpolation  and  also  increases  the 
computation  time.  Also,  density  of  the  polar  data  becomes  sparser  as  one  gets 
farther  away  from  the  center,  so  interpolation  error  becomes  greater  for  the  higher 
frequency  components  which  results  in  some  image  degradation. 

An  alternative  algorithm  which  is  described  in  the  following  sections  requires 
no  interpolation  and  uses  1-D  FFT  and  an  additional  algorithm  making  use  of  the 
well  known  “Fourier  Slice  Theorem”  (or  “Projection  Slice  Theorem”)  and  back- 
projection  techniques  [8,9,10,3]  can  be  used  to  improve  the  quality  of  the  images. 


2.2.1  Fourier  Slice  Theorem 


Let  h(x,y)  be  a  2-D  function  and  assume  its  2-D  Fourier  transform  exists. 

(2.10) 


M*,y)  & 


1 1  and  r  are  the  rotated  coordinates  as  defined  in  Equations  (2.2a)  and  (2.2b).  Let 
pe{u)  be  the  projection  of  h(u,v)  at  angle  0,  defined  by 


Pe(u)  =  f  P{x,y)dv, 

J  —  oo 

=  /  P(  u  cos  9  —  v  sin  0,  u  sin  0  +  v  cos  9)dv, 

J  —  OO 

and  assume  its  Fourier  transform,  Pg(^)  exists. 

P„{lo)  =  FT{pe(u)}  =  r  Pe(u)e~iu,u du 

J  —  OO 

/oo  roc 

/  h(u  cos  9  —  v  sin  9,  u  sin  9  +  v  cos  9)e~3uu dudv. 

-oo  J  —  oo 

Transforming  back  to  the  unrotated  coordinate  system, 

Pe=  r  f°°  h(x,y)e~Mxco*°+ySm9)dxdy,=  ff(u>  cos  fi,  cos  0). 

J  —  OO  J  -  OO 


(2.i  i ; 


(2.12) 


(2.13) 


or 


Pe{f)  =  H(fcos9Jsin9).  (2.14) 

Hence,  if  the  (u,v)  pair  is  rotated  by  an  angle  9 ,  the  Fourier  transform  of  the 
projection  is  equal  to  the  2-D  Fourier  transform  of  the  object  along  a  line  rotated 
by  9.  This  result  leads  to  the  “Fourier  Slice  Theorem”  which  is  stated  as  [9,10]:  “the 
Fourier  transform  of  a  parallel  projection  of  an  image  h(x,y )  taken  at  an  angle  9 
gives  a  slice  of  the  2-D  transform  )  subtending  an  angle  9  with  the  u  axis." 

In  other  words,  the  Fourier  transform  of  Pg(  u)  gives  the  values  of  h(fx,fy)  along 
the  line  /  in  Figure  2.2. 

The  above  theorem  indicates  that  by  taking  projections  of  an  image  at  angles 
9j(j  =  l,2,...,n)  and  Fourier  transforming  each  of  these  projections,  values  of 
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Figure  2.2:  Projection  of  a  2-D  function  and  the  slice  of  its  Fourier  transform. 

^(/rt/y)  can  be  determined  along  the  radial  lines  shown  in  Figure  2.2.  If  H[fx>fy) 
is  known  at  a  sufficiently  large  number  of  points  on  the  (fx,fy)  plane  by  taking  a 
sufficiently  large  number  of  projections,  then  the  2-D  image  of  the  object  can  be 
reconstructed  by  inverse  Fourier  transform  techniques. 

2.2.2  Image  Reconstruction 

Recalling  the  formula  for  inverse  Fourier  transform, 

h(x,y)  =  r  r  H(fxJy)e^I+f^dfxdfy.  (2.15) 

J  -  oo  J  -oo 

Defining  a  polar  coordinate  system  (/,  #)  by 

fz  =  f  cos  0,  (2.16a) 


fy  =  /  sin  9 , 


(2.16b) 


du  d V  =  fdfde 


(2.17) 


sin(0  +  7r)  =  —  sin#  (2.21b) 

the  integral  of  Equation  (2.19)  can  be  rewritten  as, 

h(x,y)  =  H ( f ,  9)\f\e}2"^T  cos$+v df  dO 

=  j’  [/°°  P,\/n\de.  (2.22) 

where  t  =  x  cos  9  +  y  sin  9. 

The  inner  integral  represents  a  one-dimensional  inverse  Fourier  transform  of  the 
product  of  the  projection  p$(f)  and  |/|.  It  thus  corresponds  to  a  filtered  function. 
Summing  the  projections  for  each  aspect  angle,  9 ,  an  estimate  for  h(x,y)  can  be 


obtained.  Substituting  y(.r,  y)  for  h(x,y)  in  Equation  (2.22),  image  reconstruction 
procedure  can  be  summarized  as  follows:  to  every  point  (x,y)  on  the  image  plane, 
there  corresponds  a  value  of  t  —  x  cos  0  4-  y  sin#  for  a  given  value  of  0  and  the 
filtered  projection  contributes  to  the  reconstruction  with  its  value  of  1.  Making  use 
of  Equation  (2.22),  the  new  algorithm  substitutes  for  the  2-D  Fourier  transform. 
Since  it  consists  of  only  a  1-D  Fourier  transform  plus  an  addition  algorithm,  it 
improves  the  quality  of  the  images  by  avoiding  the  errors  due  to  interpolation  of 


the  data  points. 


2.3  Limitations 


From  a  practical  point  of  view,  it  is  not  possible  to  have  either  a  continuous 
frequency  spectra  or  infinite  aspect  angles  around  the  target.  Therefore,  the  quality 
of  the  images  produced  also  depends  on  the  sampling  rates  of  both  frequency  and 
aspect.  On  the  other  hand,  in  order  to  represent  original  continuous  waveforms 
sampling  rates  must  satisfy  some  criterion  which  will  be  discussed  in  the  following 


sections. 


2.3.1  Frequency  Sampling  Criteria 


Discrete  data  sampling  at  the  Nvquist  rate  for  reconstruction  of  continuous 
waveforms  is  well  known.  From  the  sampling  theorem  [11,12]:  “A  time  limited 
signal  can  be  reproduced  from  its  discrete  values,  if  it  is  sampled  over  the  complete 
frequency  domain  using  the  Nvquist  rate."  Where,  /,,  the  sampling  rate  should 
satisfy  f,  <  1/27"  and  signal  is  limited  to  time  ±T. 


Then, 


w.  t  rr:y'w{''Tr  n”. 

L-~’  \  7  /  icT  n  n 


where  samples  are  taken  at  ic  ’y  .  oc.  ■  n  <  oc .  To  employ  the  sampling 


A  ii  I  l/l  ■  I  Mi  1  A  A  ft  <  Aj  1 ^ a1,  1.  >  .  a.  i  j'J ■* 


1  *>*  ■  *  ■•»*>» 
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theorem,  the  signal  must  be  limited  either  in  time  or  frequency.  In  our  cases, 
targets  are  limited  in  space,  or  equivalently  signals  are  limited  in  time.  On  the 
other  hand,  measured  spectra  is  also  frequency  limited.  Therefore,  target  images 
are  expected  to  repeat  in  space. 

Applications  of  sampling  theorem  for  general  N-dimensional  cases  has  been  stud¬ 
ied  by  Fok  [13]  and  the  proposed  frequency  increment  for  the  2-D  case  is 


A  /< 


(2.24) 


J~2(l  +  k)D'  K  ’ 

where,  D  is  the  maximum  dimension  of  the  object,  C  is  the  speed  of  light,  and  K 


is  some  safety  factor. 


2.3.2  Number  of  Projections  Required 


Given  that  one  can  only  obtain  a  finite  number  of  projections,  Ar,  what  should 
the  minimum  value  of  N  be?  This  depends  to  a  great  extent  upon  what  is  known 
a  priori  about  the  signal  being  reconstructed  and  how  much  detail  we  expect  to 
learn  about  it.  For  example,  if  it  is  known  a  priori  that  an  unknown  signal  is 
circularly  symmetric  (like  scattering  from  a  sphere),  all  its  projections  are  identical 
and  only  one  of  them  is  necessary  (and  will  be  sufficient)  to  reconstruct  g(x,y).  On 
the  other  hand,  if  absolutely  nothing  is  known  beforehand  about  g(x,y)  then  an 
infinite  number  of  projections  will  be  required.  To  determine  the  minimum  number 
of  projections  required,  the  angular  increment  which  satisfies  the  Nyquist  criterion 
is  given  by  Mensa  et  al.,  [5] 
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3.  APPLICATIONS 


Up  to  this  point,  it  has  been  shown  how  a  2-D  image  of  a  scattering  target  can 
be  reconstructed  in  the  spatial  coordinates.  On  the  other  hand,  scattering  charac¬ 
teristics  of  most  targets  depend  upon  the  polarization  of  the  incident  wave  which 
requires  one  more  dimension  to  be  introduced,  namely  polarization,  to  make  the 
images  more  identifiable.  Early  studies  on  polarization  properties  of  the  scattered 
waves  have  been  carried  on  by  Sinclair  [14],  Kennaugh  [15]  and  Graves  [16].  More 
recently,  this  property  of  scattered  waves  has  been  studied  by  Boerner  et  al.  [17], 
and  Kostinski  et  al.  [18].  Also  in  most  of  the  papers  concentrated  on  the  scattering 
from  different  geometries  using  different  approximation  techniques,  such  as  GTD, 
one  can  find  the  formulation  of  the  scattered  field  for  that  specific  geometry  for 
different  polarizations. 

In  this  report,  images  of  some  geometrically  simple  (sphere,  tilted  square  plate, 
tilted  ellipsoid)  and  more  complicated  (an  aircraft)  targets  using  either  measured  or 
simulated  data  were  illustrated  and  analyzed.  Section  3.1  briefly  discusses  measure¬ 
ment,  approximations  used  to  generate  the  simulated  data  and  the  data  processing 
technique  used  for  imaging.  Section  3.2  illustrates  the  images  produced  including 
the  discussions  related  to  them. 

3.1  Data  Processing 

Scattering  measurements  of  several  targets  were  made  using  The  Ohio  State 
University  Compact  Radar  Cross  Section  Measurement  Range  [6].  Measured  raw 
data  are  calculated  using  the  calibration  equations  given  by  the  same  reference. 
After  calibration  and  windowing  processes  the  measured  data  become  ready  for  the 
imaging  process.  Simulated  data  were  generated  using  approximations  based  on 
uniform  asymptotic  techniques  formulated  in  terms  of  the  Geometrical  Theory  of 
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Diffraction  (GTD)  [7] .  The  imaging  process  can  be  summarized  as  follows: 

1.  Multiply  the  frequency  domain  data  by  the  frequency  (Equation  (2.21)). 

2.  Inverse  Fourier  transform  the  filtered  frequency  data  of  Step  1. 

3.  Determine  the  image  size  which  gives  the  size  of  the  required  time  domain 
window. 

4.  Average  a  certain  number  of  points  on  the  time  domain  data  (this  number  is 
determined  by  the  time  domain  window  of  step  3)  so  that  the  image  can  be 
represented  by  a  n  x  w  matrix  (100  x  100  is  used). 

5.  Calculate  the  discrete  coordinates  on  the  image  plane  (a  100  x  100  matrix). 

6.  For  each  discrete  point  on  the  image  plane,  sum  the  contribution  of  each  time 
domain  response  by  rotating  the  time  domain  signals  by  9. 

7.  If  more  than  one  polarization  is  required  at  the  same  time,  go  back  to  Step  5 
for  the  new  polarization. 

A  frequency  range  of  2-18  GHz  has  been  used  in  most  of  the  measurements. 
For  square  plate  applications,  two  different  scale  models  (ratio  <  8  :  1)  have  been 
combined  to  widen  the  bandwidth  to  0.01  to  18  GHz  such  that  low  frequency  effects 
can  also  be  observed.  Simulated  data  spans  the  same  frequency  range  in  order  to 
directly  compare  measured  and  calculated  images. 

Two  linear  polarizations,  vertical  (VP)  and  horizontal,  are  used  to  demonstrate 
the  polarimetric  scattering  characteristics  of  the  targets  of  interest.  A  Tekronix  429 
color  CRT  display  and  hard  copy  unit  are  used  to  display  the  images.  Two  different 
sets  of  colors  are  assigned  to  encode  the  amplitude  for  each  polarization.  The 
computer  programs  CIMAG2  [22]  and  CLRPL  (Appendix  I)  are  used  to  generate 
the  images  and  to  plot  them  on  the  color  display,  respectively. 
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3.2  Results 


The  first  target  of  interest  is  a  sphere  whose  impulse  response  and  images  with 
different  spatial  windows  are  shown  in  Figures  3.3  through  3.5.  Because  it  is  a 
geometrically  simple  object,  its  scattering  characteristics  are  well  known  [21,19], 
and  also  it  is  aspect  and  polarization  independent.  Hence,  it  is  a  useful  introductory 
step  to  start  imaging  studies  with  a  sphere  to  examine  the  imaging  technique  used 
and  understand  different  features  on  the  images.  Figure  3.4  shows  the  image  of 
the  sphere  with  a  360°  spatial  window  with  1°  increments.  The  image  of  the  same 
sphere  with  a  90°  spatial  window  is  shown  in  Figure  3.5.  The  frequency  data  used  to 
generate  these  images  range  from  2  to  18  GHz  in  10  MHz  steps.  Examining  Figure 
3.4,  the  red  colored  ring  with  the  highest  amplitude  (or  the  red  colored  arc  in  Figure 
3.5)  corresponds  to  the  outer  radius  of  the  sphere  which  is  related  to  the  specular 
response.  Green  colored  rings  outside  the  sphere  with  relatively  lower  amplitudes 
correspond  to  the  higher  order  scattering  mechanisms  like  creeping  waves  around 
the  sphere.  Three  concentric  circles  inside  the  red  one  aTe  related  to  the  ripples 
around  the  specular  in  Figure  3.3.  The  same  mechanisms  discussed  above  can  be 
observed  in  Figure  3.5.  Since  only  a  limited  spatial  window  is  used,  the  figure  is 
not  symmetric  anymore.  Therefore  the  delayed  response  related  with  the  creeping 
waves  which  is  represented  by  the  green  arc  facing  the  red  one  can  be  more  clearly 
identified  in  this  figure. 

Figures  3.7  through  3.14  belong  to  a  1  ft.  tilted  (45°  from  horizontal)  square  plate 
using  either  measured  or  simulated  data  with  different  spatial  windows  (90°  and 
360°)  and  linear  polarizations  (scattering  mechanisms  related  to  flat  plate  structures 
has  been  analyzed  and  discussed  in  detail  by  Sitka  [20]).  Frequency  data  spans  a 
range  of  0.01  -  18  GHz  with  10  MHz  steps,  which  was  achieved  by  combining  the 
frequency  data  belonging  to  2  different  scale  modeled  plates  for  the  measured  data. 
An  impulse  response  waveform  belonging  to  one  of  these  combined  data  is  shown 
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in  Figure  3.6.  First,  comparing  Figures  3.9  and  3.13,  which  are  the  images  of 
the  same  plate  with  a  spatial  window  of  90°  using  different  linear  polarizations, 
a  stronger  response  is  observed  from  the  back  edge  in  VP  while  the  front  edge 
is  more  dominant  in  HP,  because  the  surface  waves  cannot  propagate  along  the 
plate  and  contribute  the  scattering  at  the  back  edge  in  the  HP  case.  A  second 
order  mechanism  corresponding  to  the  waves  traveling  along  the  plate  is  seen  as  a 
delayed  response  in  Figure  3.9  (VP)  which  is  coded  by  a  yellow  green  line  behind 
the  back  edge  and  it  does  not  show  up  in  Figure  3.13  because  of  the  same  reason 
discussed  above.  Response  around  90°  incidence  (side  edge  on  the  right)  is  almost 
identical  in  both  figures  as  expected.  Strong  responses  at  the  corners  are  observed 
on  both  polarizations.  Now,  comparing  Figures  3.9  and  3.10  or  3.13  and  3.14  a 
good  agreement  is  observed  between  the  measured  and  calculated  images  except 
that  the  second  order  mechanism  observed  on  the  measured  image  isn’t  seen  in 
the  calculated  ones  because  simulated  data  are  generated  by  a  first  order  UTD 
approximation.  The  same  discussions  hold  for  Figures  3.7,  3.8,  3.11  and  3.12,  the 
images  with  full  360°  spatial  window.  One  advantage  of  using  a  90°  spatial  window 
is,  that  figures  are  not  symmetric  any  more  and  it  is  easier  to  identify  the  delayed 
responses  (it  also  requires  far  less  data).  In  all  of  the  square  plate  images  discussed 
above,  the  increment  for  the  aspect  angle  (A0)  used  is  15°  which  does  not  satisfy  the 
sampling  criterion  and  causes  some  image  degradation.  Also,  while  combining  some 
of  the  frequency  data  belonging  to  different  scaled  modeled  plates,  some  alignment 
problems  caused  the  combined  data  to  be  slightly  blurred  which  produced  some 
noise  in  the  measured  images  too. 

Figures  3.14  -  3.20  belong  to  a  tilted  ellipsoid  with  different  elevations  and  two 
different  linear  polarizations.  Frequency  data  spans  a  1-12  GHz  range  with  10  MHz 
steps.  The  increment  for  the  aspect  angle  is  30°  which  doesn’t  satisfy  the  sampling 
criteria.  When  Figures  3.14  -  3.16  or  3.1  7  -  3.20  are  examined,  the  size  of  the  figures 
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gets  smaller  as  the  tilt  angle  becomes  larger  because  the  projected  area  becomes 
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Figure  3.4:  Impulse  response  of  a  1  ft.  tilted  square  plate.  Elevation  —  45°;  aspect 
=  90°;  frequency  =  0.01  -  18  GHz.  Result  of  combination  of  two  differently  scaled 
modeled  plates  (ratio  8:1)  with  2-18  GHz,  measured  data. 
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Figure  3.12:  Image  of  a  1  ft.  tilted  square  plate  with  a  90°  spatial  window. 
Elevation  —  45°;  IIP;  calculated  data. 
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Figure  3.13:  Image  of  an  ellipsoid.  Elevation  0°;  VP. 
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Figure  '{.21:  Image  of  a  Boring  707  aircraft,  IIP. 
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4.  CONCLUSIONS 
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It  was  shown  that  the  2-D  image  of  the  target  is  related  to  its  2-D  frequency 
spectra  by  a  Fourier  transform  relationship.  However,  the  frequency  domain  data 
is  sampled  on  a  polar  raster  and  needs  to  be  transformed  into  a  rectangxdar  one. 
The  imaging  algorithm  described  in  this  report  substitutes  for  conventional  2-D 
Fourier  transform  and  avoids  the  errors  due  to  interpolation  by  using  a  1-D  inverse 
Fourier  transform  and  a  summation  procedure.  Hence,  it  improves  the  quality  of 
the  images. 

Images  are  polarization  dependent  because  the  scattering  properties  of  most 
scatterers  are  polarization  dependent.  Use  of  polarization  information  as  a  new 
dimension  in  images  makes  them  more  uniquely  identifiable.  Use  of  color  helps  dif¬ 
ferent  scattering  centers  related  to  different  polarizations  to  be  identified,  especially 
when  the  images  corresponding  to  different  polarizations  are  combined  in  the  same 
frame  which  may  be  a  part  of  the  future  research  as  a  continuation  of  the  work 
described  here. 

A  sphere,  tilted  square  plate,  tilted  ellipsoid  and  a  Boeing  707  aircraft  are  the 
targets  investigated  in  this  research.  The  geometrically  simple  ones  help  to  under¬ 
stand  the  basic  scattering  characteristics.  The  square  plate  was  studied  both  with 
measured  and  theoretical  data.  Results  are  found  to  be  in  good  agreement,  and 
polarization  is  found  to  be  a  useful  tool  for  target  recognition  and  identification 
with  imaging.  More  work  is  required  to  be  done  to  improve  the  aircraft  images. 
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THE  COMPUTER  PROGRAM  CLRPL 


THE  SOFTWARE  CODE  TO  DISPLAY  IMAGES  ON  THE  TEKTRONIX  429 
COLOR  CRT  DISPLAY. 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

cc  cc 

CC  CLRPL  CC 

cc  cc 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

c 

C  THIS  PROGRAM  DEFINES  A  PLOTTING  SOFTWARE  CODE  TO  PLOT 

C  IMAGE  OF  AN  OBJECT  WITH  POLARIZATION  INFORMATION 

C 

C  IMAGE  DATA  SHOULD  BE  CODED  BY  A  100*100  MATRIX  AND 

C  IT  SHOULD  BE  CONTAINED  IN  THE  MATRIX  ’ ARRAY ( 1 00 , 1 0 0 ) * 

DIMENSION  ICOL( 100, 100 ) , I COLA ( 100 ) 

DIMENSION  GX( 5) ,GY<  5) ,CBX( 5 ) ,CBY( 5 ) , XP ( 5 ) , YP ( 5 ) , XT (  1 0 ) , YT ( 1 0 ) 
DIMENSION  ARRAY( 100 , 100 ) , COORD( 5 ) 

INTEGER  COLI , WXSTI D 
REAL*4  L 

LOGICAL  FINISHED 
CHARACTER* 4 0  FILNM 
CHARACTER *7  NUMB 

C  THE  GX  MATRIX  IS  THE  X  COORDINATES  OF  THE  COLOR  GRID 

DATA  GX/. 30, .85, .85, .30, .30/ 

C  THE  GY  MATRIX  IS  THE  Y  COORDINATES  Or  THE  COLOR  GRID 

DATA  GY/. 12, .12, .67, .67, .12/ 

C  THE  CBX  MATRIX  IS  THE  X  COORDINATES  OF  THE  COLOR  BOX 

DATA  CBX/. 30, .85, .85, .30, .30/ 

C  THE  CBY  MATRIX  IS  THE  Y  COORDINATES  OF  THE  COLOR  BOX 

DATA  CBY/. 70, .70, .72, .72, .70/ 

C  THE  XT  MATRIX  IS  THE  X  COORDINATES  OF  THE  TIK  MARKS 

DATA  XT/. 3, .44, .57, .71, .85, .26, .26, .26, .26, .26/ 

C  THE  YT  MATRIX  IS  THE  Y  COORDINATES  OF  THE  TIK  HARKS  ON  THE  Y  AXIS 

DATA  YT/.08, .08, .08, .08, .08, .12, .26, .40, .54, .67/ 


PI-4 . *ATAN ( 1 . ) 

KWK-412505 

WKSTID-1 

C 

C  YMTCL — Y  Value  of  Magnitude  Top  Center  Line 

YMTCL— CBY ( 3 )  +  . 0  2 
C 

C  ENTER  THE  ARRAY  TO  BE  PLOTTED 

WRITE(6, «) 'ENTER  FILE  NAME' 

READ (5,100)  FILNM 
100  FORMAT (A40 ) 

C 

C 

C  POLARIZATION  DEPENDENT  PLOT? 

WR I TE ( 6 , * )  'DO  YOU  NEED  TO  USE  POLARIZATION  AS  A  PARAMETER?’ 
WR I TE ( 6 , * ) ’ ( Y-l  N-0) ’ 

READ ( 5 , * )  POL 
ir(POL.EQ.l)  THEN 
C  ENTER  THE  POLARIZATION 

WRITE(6, *) ’ENTER  POLARIZATION  (VP-1, HP-2)’ 

READ! 5,*)  PL 
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c 


ELSE 
END  IF 


OPEN ( UNIT-1 , FI LE-FILNM, STATUS- 'OLD' , FORM- ' UNFORMATTED' ) 
C 

C  ENTER  THE  IMAGE  SIZE  AND  PERIOD  OF  THE  TIME  SIGNAL 

READ ( 1 )  I  MGS  Z , PER 
C 

C  ENTER  THE  IMAGE  ARRAY 

DO  1000  1-1,100 
DO  1000  J-1,100 
READ( 1 ) ARRAY ( I , J ) 

ARRAY ( I , J ) -ABS ( ARRAY ( I , J ) ) 

1000  CONTINUE 

C 

C  NORMALIZE  THE  ARRAY  VALUES 

C 


CALL  SEARCH ( ARRAY , AMAX , AMIN ) 

WRITE( 6,* ) 'MAX-' , AMAX , ' MIN- AMIN 

WRITE( 6, * ) ' ENTER  THE  DESIRED  MAX . , AND , MIN . ' 

READ ( 5 , * )  ANMAX , ANMIN 
DO  20  1-1,100 
DO  20  J-1,100 

IF  (ARRAY) I , J) .GT. ANMAX)  ARRAY ( I , J ) -ANMAX 
IF  ( ARRAY ( I , J ) . LT . ANMI N )  ARRAY) I , J ) -ANMIN 
ARRAY ( I , J ) - ( ARRAY ( I , J ) -ANMIN )/ ( ANMAX-ANMIN ) 

20  CONTINUE 

C 

C  OPEN  GKS  ERROR  FILE 

CALL  GOPKS ( 6 , S000 ) 

C  FIND  CONNECT  ID 

999  CALL  GKHGCI ( ' ESL  4129 ' , JERROR , RCONID ) 

I F ( J ERROR . NE . 0 )  THEN 

WRITE<6,*)  'Can  not  be  a  connection  ID' 
wri TE ( 6 , *  )  'Would  you  like  to  wait  ?  (Y-l)' 

READ) 5, * )ANS 
IF(ANS.NE.l)  THEN 
STOP 
ELSE 

WRITE ( 6 , * )  'Enter  1  when  ready' 

READ) 5,*)  ANS 
GO  TO  999 
END  IF 
END  IF 

C  OPEN  WORKSTATON  #1 

CALL  GOPWR) 1 .RCONID, EWE) 

C  ACTIVATE  WORRSTATION  fl 

CALL  GACWK(l) 

C  SET  THE  WORRSTATION  WINDOW/VIEWPORT-FULLSCREEN 

C  GET  MAX  X  AND  Y 

RUNITS-0 

CALL  GQDSP) RWR , RERROR , RUNITS , XS I ZE , YSI ZE , RRASX , RRASY ) 
CALL  GSWRWN (1,0. , 1 , , 0 . , YSI ZE/XSIZE ) 

CALL  GSWRVP) 1, 0  ...  343,0. , .274 ) 


GENERATE  THE  COLOR  INDICES 


IF(POL.EQ.O)  THEN 
DO  3  COLI-1 , 100 
L-.5 
S-l 
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IT(COLI .LE.5)L-0 

H-95+2.6*COLI 

CALL  HLSRGB(H,L,S,R,G,B) 

CALL  GSCR(1,1+C0LI,R,G,B) 

ICOLA( COLI j-l+COLI 

3  CONTINUE 

ELSE 

DO  4  COLI-1,100 
L-.5 
S-l 

FC-25*(4.8)»*(COLI/100. ) 

I F ( COLI . LE . 5 )  L-0 
IF(PL.EQ.l)  H-225-FC 
IF(PL.EQ.2)H-U0-FC 

I F ( PL . EQ . 2 . AND . COLI . LT . 2  5 . AND . COLI . GT . 5 ) L- . 7  5 

IF(H.LT.O)  H-H+360 

CALL  HLSRGB(H,L,S,R,G,B) 

CALL  GSCR(l,l+COLI,R,G,B) 

ICOLA( COLI )-l+COLI 

4  CONTINUE 
END  IF 

GENERATE  COLOR  CODE  FOR  X  AND  Y  COORDINATED  RCS  LEVELS 
DO  30  1-1,100 
DO  30  J-1,100 

I COL (  I ,  J  ) -ARRAY ( I ,  J  )  *99  +  2 
0  CONTINUE 

PLOT  COLOR  LABEL  USING  CELL  ARRAY 
CALL  GCRSG(SEG) 

CALL  GCA ( CBX ( 1 ) ,CBY( 1 ) ,CBX( 3 ) ,CBY<  3), 100, 1,1, 1,100,1, I COLA) 
CALL  GSPLCI(l) 

CALL  GPL ( 5 , CBX , CBY ) 


LABEL  COLOR  BAR 
CALL  GSTXAL(2,3) 

CALL  GSTXP(O) 

CALL  GSCHXP (1.25) 

CALL  GSCHSP ( 1 . ) 

CALL  GTX(GX( 1 ) , YMTCL,  '  0.0'  ) 

CALL  GTX( .44, YMTCL, ’0.25' ) 

CALL  GTX( . 57 , YMTCL, '0.50' ) 

CALL  GTX( .71, YMTCL, '0.75' ) 

CALL  GTX(GX( 2 ) , YMTCL, '1 .0’ ) 

CALL  GSCHXP ( 1 . ) 

CALL  GSCHSP(  1 .  ) 

PLOT  COLOR  MATIX  USING  CELL  ARRAY 
TEST-10 

IF  ( TEST . EQ . 0 ) GO  TO  21 

CALL  GCA(GX(1),GY(1),GX(3),GY(3),100,100,1,1,100,100,ICOL) 
CALL  GSPLCI(l) 

CALL  GPL( 5 , GX , GY ) 

DRAW  GRID  LINE  ADJACENT  TO  EACH  AXES 

XP( 1 )-. 30 
XP( 2 )- . 85 
YP ( 1 ) • . 08 
YP( 2 )- . 08 
CALL  GPL(2,XP,YP) 

XP< 1 )-. 26 
XP( 2 )- . 26 
YP ( 1 )  -  .  1 2 
YP ( 2  )  -  .  67 
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CALL  GPL ( 2 , XP , YP ) 

C  PRINT  THE  TIK  KARRS 

CALL  GSMK ( 2 ) 

CALL  GSMKSC(l) 

CALL  GSPMCI(l) 

CALL  GPM (10, XT , YT ) 

C 

C  SHOW  THE  ASPECT  ANGLES 

XP( 1 ) - . 57 
XP( 2 )-. 57 
YP ( 1 )-  .  16 
YP( 2 )-. 19 
CALL  GPL ( 2 , XP , VP) 

XP< 1 )-.S65 

XP ( 2  )  - . 57 

XP ( 3 ) - . 575 

YP ( 1 )  —  - 1 8  5 

YP ( 2 ) - - 1 9 

YP<  3)-. 185 

CALL  GPL ( 3 , XP , YP ) 

XP( 1  )-.80 
XP ( 2  )• .  77 
YP ( 1  )  - .  40 
YP( 2  )  - . 40 
CALL  GPL ( 2 , XP , YP ) 

XP( 1 )- . 775 
XP ( 2 ) - . 77 
XP ( 3 ) - . 775 
YP ( 1 )-. 405 

YP( 2 )-. 4 

YP( 3 )- . 395 
CALL  GPL ( 3 , XP , YP ) 

CALL  GTX( .57, .14, '0  DEG.’) 

CALL  GTX( .57, .02, 'TIME  IN  NANOSECONDS') 

CALL  GSTXAL(2,1J 
CALL  GSTXP(O) 

CALL  GSCHUP ( 1 . ,0. ) 

CALL  GTX( .82, .4, '90  DEG.') 

CALL  GTX( .20, .40, 'TIME  IN  NANOSECONDS’) 

C 

C  riGURE  OUT  THE  COORDINATES  ON  EACH  AXES 

CALL  MARKS ( IMGSZ, PER, COORD) 

C 

C  PRINT  THE  COORDINATES 

X-.30 

CALL  GSCHUP( 0 . ,1. ) 

DO  40  1-1,5 

WRITE ( NUMB, TMT-' ( P6 . 2 ) ’ )COORD( I ) 

CALL  GTX ( X , .06, NUMB) 

40  X-X-f.14 

CALL  GSCHUPU.  ,0.  ) 

Y-.12 

DO  50  1-1,5 

WRITE ( NUMB , FMT- ’ (T5. 2 ) ’ ) COORD ( I ) 

CALL  GTX( . 23, Y, NUMB) 

50  Y-Y+.14 

C 

C 

C  COLOR  1  -  WHITE,  0  -  BLACK  (ON  PLOTTER,  REVERSED  ON  SCREEN) 

C 

21  CALL  GCLSG(SEG) 

C 

WRITE( 6 , * ) 'Enter  return  to  finish . ' 

READ (5,1 ) FINISHED 


46 


sS&S 


y  s 


1  FORMAT! Al) 

CALL  GCLRWK ( WKSTI D ,  1 )  (Clears  the  screen  on  TEXTRONIX 

C  DEACTIVATE  THE  WORK  STATION 

CALL  GDAWK(l) 

C 

C  CLOSE  THE  WORK  STATION 

CALL  GCLWK(l) 

C  CLOSE  THE  SYSTEM 

CALL  GCLKS 
STOP 
END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CC  THIS  ROUTINE  MAKES  THE  TRANSFORMATION  BETWEEN  COLOR  SYSTEMS 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  HLSRGB(H,L,S,R,G,B) 

REAL  H,L,S,R,G,B,Ml,M2 
IF  (L  .LE.  .5)  THEN 
M2-L* ( 1+S ) 

ELSE 

M2-L+S-L*S 
8  END  IF 

Ml-2*L-M2 

B-rgb_value(Ml ,M2,H+120) 

R«  r gb_value ( Ml , M2  ,  H ) 

G-rgb_value(Ml ,M2 .H-120 ) 

RETURN 

END 

FUNCTION  rgb_value ( Nl , N2 , HUE ) 

REAL  rgb_value , Nl , N2 , HUE 
IF  (HUE  .GT.  360)  THEN 
HUE-HUE-360 
END  IF 

IF  (HUE  .LT.  0)  THEN 
HUE-HUE+360 
END  IF 

IF  (HUE  .LT.  60)  THEN 

rgb_value-Nl+ ( N2-N1 ) *HUE/60 
ELSE  IF  (HUE  .LT.  180)  THEN 
rgb_value-N2 

ELSE  IF  (HUE  .LT.  240)  THEN 

rgb_value-Nl-MN2-Nl  )*( 24 0-HUE )/60 
ELSE 

rgh_va lue-Nl 
END  IF 
RETURN 
END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  THIS  ROUTINE  FINDS  THE  MAX.,  AND  MIN.  OF  A  100*100  ARRAY 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C 

SUBROUTINE  SEARCH ( ARRAY , AMAX , AMIN ) 

DIMENSION  ARRAY (100,100) 

AMAX--1000 
AMIN-1000 
DO  1  1-1,100 

DO  1  J-1,100 
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I F ( ARRAY ( I , J) . GT.AHAX)AHAX-ARRAY ( I , J) 

I F ( ARRAY (  I ,  J  ) . LT . AM I N ) AMI N« ARRAY ( I ,3) 

1  CONTINUE 

RETURN 
END 

ccccccccccccccccccccccccccccccccccccccccecccccccccccceccccccccccccccccccccccccc 

C  THIS  ROUTINE  CALCULATES  THE  DIVISIONS  ON  THE  AXES 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


SUBROUTINE  HARKS { I HGSZ , PER , COORD ) 
D I  HENS  I ON  COORD(S) 

SIZE-PER* IHGSZ/4096 . 

COORD ( 1 ) «S I ZE/-2 . 

COORD ( 2)-SIZE/-4. 

COORD ( 3 )-0 
COORD ( 4 )-SIZE/4 . 

COORD( 5)«SIZE/2. 

RETURN 

END 
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